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Abstract 

We analyze the global phase diagram of a Maier-Saupe lattice 
model with the inclusion of disorder degrees of freedom to mimic a 
mixture of oblate and prolate molecules (discs and cylinders). In the 
neighborhood of a Landau multicritical point, solutions of the statis- 
tical problem can be written as a Landau-de Gennes expansion for 
the free energy. If the disorder degrees of freedom are quenched, we 
confirm the existence of a biaxial nematic strucure. If orientational 
and disorder degrees of freedom are allowed to thermalize, this bi- 
axial solution becomes thermodynamically unstable. Also, we use a 
two-temperature formalism to mimic the presence of two distinct re- 
laxation times, and show that a slight departure from complete ther- 
malization is enough to stabilize a biaxial nematic phase. 

1 Introduction 

The transition between a uniaxial nematic structure and an orientation- 
ally disordered phase is perhaps the most investigated and best character- 
ized phase transition phenomenon in liquid crystalline systems [Ij. This 
weak first-order transition is quite well described by the mean-field theory 
of Maier and Saupe [2][3][1], which can also be formulated in terms of a 
fully-connected statistical lattice Hamiltonian [5]. The existence of a biaxial 
nematic nematic phase, however, and the transitions between different types 



of nematic structures, which have been proposed on the basis of phenomeno- 
logical calculations for systems with intrinsically biaxial molecular groups [B], 
turned out to be much more difficult to characterize experimentally [8]. Al- 
though there have been some recent reports of a biaxial nematic structure in 
thermotropic liquid crystalline systems formed by bent-core or boomerang- 
shaped molecules ^ , a biaxial phase has been first characterized in the phase 
diagram of a lyotropic liquid mixture [IHIIII], which is better represented by 
a lattice model of shape-disordered uniaxial molecules. We then revisit the 
problem of a Maier-Saupe lattice model, with the inclusion of extra degrees 
of freedom to mimic a mixture of oblate and prolate molecules (discs and 
cylinders). Slightly different versions of this problem have been analyzed by 
different authors , with some conflicting results [121 [131 El US] • According 
to the conclusions of a recent calculation for a shape-disordered Maier-Saupe 
model with restricted orientations [16], we point out that the introduction 
of two sets of degrees of freedom opens the possibihty of choosing different 
relaxation times, with different outcomes for the thermodynamic stability of 
a biaxial nematic structure. 

We then formulate and analyze a Maier-Saupe lattice model for a mixture 
of discs and cylinders. This problem includes orientational (quadrupolar) and 
shape-disordered (discs and cylinders) degrees of freedom, which might be 
associated with different relaxation times. Solutions can be obtained by the 
application of well-known methods of statistical mechanics. First, we treat 
the case of flxed (frozen) disorder, as in a typical problem of a disordered 
solid state system. In this quenched case, disorder (shape) degrees of free- 
dom are flxed, frozen, while the orientational degrees of freedom are allowed 
to thermalize during experimental times. We then treat the case of annealed 
(thermalized) disorder, in which both orientational and disorder degrees of 
freedom are allowed to reach thermodynamic equilibrium during experimen- 
tal times, and which is certainly more adequate to explain the behavior of 
a liquid system. In agreement with previous calculations for similar models, 
in the quenched case we show the existence of a biaxial nematic structure 
for typical distributions of shape-disorder [TB]. In the thermalized case, how- 
ever, there is a biaxial nematic solution of the model equations, but it turns 
out to be thermodynamically unstable (in agreement with older calculations 
by Palffy-Murhoray and collaborators [12] )• We then resort to a formalism 
based on two heat reservoirs, at distinct temperatures, which is conceived 
to describe an intermediate situation, between fully flxed and fully thermal- 
ized distributions of disorder variables [T7]. As in the calculations for the 
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Maier-Saupe model with restricted orientations [12], we show that a small 
difference of temperatures, which is equivalent to a slight departure from 
complete thermalization, is already sufficient to produce a stable biaxial ne- 
matic phase. 

It should be pointed out that we write closed-form solutions of the statis- 
tical problem, which are not restricted to the neighborhood of the transitions, 
and which can be used to draw global phase diagrams, in terms of temper- 
ature and either concentration or chemical potential. Also, it is feasible to 
extend these calculations beyond the mean-field level, as it has been done in 
the annealed case for the analogous model with restricted orientations [T8] . 
We use the model solutions to write a Landau-de Gennes expansion for the 
free energy in terms of the invariants of the tensor order parameter, but with 
model-dependent coefficients. Corresponding phenomenological expansions 
have been investigated by a number of authors |2][T9][20j, and we can use 
some asymptotic results to confirm the numerical analysis in the neighbor- 
hood of a Landau multicritical point. 



2 Disordered Maier-Saupe model 

The Maier-Saupe theory of the nematic-isotropic phase transition can be 
obtained from a statistical calculation for a fully-connected lattice model 
given by the quadrupolar Hamiltonian 

^M} = - E 4 E ^rT' (1) 

where the sum is over all pairs of lattice sites, A/N is a scaled interaction 
energy, and {S^^} is a set orientational (quadrupolar) variables, given by 

sr = \ iH< - V) , (2) 

where r^i = (n^, rv{,nf), with \ ri^i \ = 1, is a local nematic director. This mean- 
field Maier-Saupe (MS) model is known to reproduce the main features of 
the (weak) first-order transition between uniaxial nematic and disordered 
phases. 

We mimic the behavior of a binary mixture of oblate and prolate molecules 
(discs and cylinders) by introducing an additional set of (shape) degrees of 
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freedom, {Aj}, with Aj = ±1, for i = 1,...,N. Given the configurations of 
orientational and disorder (shape) degrees of freedom, {n^} and {Aj}, the 
simplest Maier-Saupe Hamihonian for this mixture of discs and cyhnders is 
given by 

l<i<j<-'V fj.,u=x,y,z 

which can also be written in the more convenient form 

2 

(4) 

where we have discarded irrelevant terms in the iV — )• oo limit. 

In a typical problem of a disordered system of solid state physics, the dis- 
ordered degrees of freedom are fixed, frozen, while the orientational degrees 
of freedom are allowed to thermalize during the experimental times [21]. In 
this fixed, quenched case, as in amorphous and glassy materials, disorder 
variables are not strictly thermodynamic. In the opposite case, which seems 
more adequate to describe liquid mixtures, both orientational and disorder 
degrees of freedom are allowed to reach thermodynamic equilibrium during 
experimental times. The fully thermalized, annealed case, is then treated 
according to the standard rules of equilibrium thermodynamics. 

In the following paragraphs, we consider quenched and annealed cases 
separately. As in the work of Henriques and Henriques [13], for a lattice 
Maier-Saupe model with restricted orientations, we confirm that there is a 
biaxial nematic structure in the quenched case. Also, we show that this bi- 
axial structure becomes thermodynamically unstable in the annealed case, 
which agrees with an older Maier-Saupe calculation for a mixture of cylin- 
ders and discs by Palffy-Muhoray and collaborators [12]. We then introduce 
the two-temperature formalism [16] [IT] in order to show that a slight depar- 
ture from complete thermalization is already sufficient to produce a stable 
biaxial nematic phase. The numerical analysis of the free energy is supple- 
mented, and confirmed, by an analysis of a Landau-de Gennes expansion in 
the neighborhood of the Landau multicritical point. 
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2.1 Quenched disorder 

Given the set of disorder variables, {K}, we write the canonical partition 
function 



Z({A.})=i;exp^ 

{51} 



fJ,,U 



N 



i=l 



(5) 



where (3 = A/ksT = 1/t is the inverse of a (dimensionless) temperature, 
/X, = 1, 2, 3 correspond to the Cartesian directions, and we are summing 
over configurations of the local (microscopic) directors. In this quenched case, 
{Aj} is a set of independent, identical, and identically distributed random 
variables, associated with a probability distribution 



N 



^({A.})=n^^(A.). 



(6) 



It is convenient to parametrize the local directors by polar coordinates, 

r^i = (sin 9i cos (pi, sin 9i sin (pi, cos 9i) , (7) 
with fli = {9i, (pi), so that 

sin^ 9i cos^ (pi — k sin^ 9i sin (pi cos (Pi sin 9i cos 9i cos 
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S^" — - I sin^ 9i sin (pi cos 
sin 9i cos 9i cos ( 



I sin^ 9i sin (pi cos ( 

sin^ 9i sin^ (pi — | sin 9i cos 6*^ sin (pi | , (8) 

cos^ 9i — ^ 



sin cos sin ( 



and the sum over orientational configurations becomes an integral over solid 
angles, 

£_ 
2N 



N 



1=1 



(9) 



In the thermodynamic limit (A'^ oo), we have the asymptotic form 

Z{{X,})r^exp[-/3Ngr,{{Xi})]. (10) 

The resulting (quenched) free energy Qq comes from an average of Qn {{Xi}) 
over the distribution of shape variables P {{Xi}), 



^ (In Z ({A,})) = ^J (U dx}j P ({A,}) In Z ({A,} , /3) , 

(11) 
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where the brackets (...) indicate disorder averages, and we are taking the 
hmit of large A^. 

The sum over the square terms in equation (|9]), can be deah with by a 
set of Gaussian identities. For example, we have 



exp 



2N 



dxii 



exp < —X 
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2N 



N 



i=l 



Xii 



(^) ^ /^?iiexp|-^Ar/3g?i + 5^/3A.^"(f^.)gii|. 



:i2) 



Taking into account the symmetry of the traceless tensor S^'^ , we introduce 
a set of six variables, gn, ^22, fe, ^12, ^is, and ^23, and write the partition 
function 



Z{{\}) 
where 



[dq] exp - E 9^ + E 1^ ^4 ' ^^^^ 

I, /i /i<!^ i=l ) 



[dq] 



I3N 



dqndq22dq33dqi2dqndq23, 



and 



(14) 
(15) 



In the thermodynamic limit, we resort to Laplace's asymptotic method, 
and invoke the law of large numbers, 

1 ^ r 

-ElnM(A„{g^J)-^(lnM(A„{g^4))= / rfAp (A) In M (A, {g^J) . 
i=i 

(16) 

We then have a self-averaged expression for the quenched free energy, 
9g = ^ [qli + 922 + + ?i2 + gi3 + to - ^ / (^) In M (A, {g^,}) , (17) 
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where the set of parameters {g^i^} come from the minima of the asymptotic 
integration, 



/ dns^^ (n) exp 



qss = / (A) d\- 



.fJ.<U 



J (in exp 



.fJ,<U 



;i8) 



for (5 = 1,2, 3, and 



/ dnS^^ (n) exp 



J dQexp 



.fJ.<U 



(19) 



for (5 < 7. The set of variables {g^i/} has a clear physical interpretation as 
the mean values of the quadrupole tensor components. In fact, if we include 
field terms in the original Hamiltonian, of the form h^^XiS^^ , with couplings 
of local quadrupoles XiS'^" to external fields h^'' , the (Gibbs) free energy will 
depend on these external fields, and the mean quadrupoles will be given by 
qf,u = -dgq/dh^^. 

Using the explicit forms of 5'^'^ (fi), given by equation ([8]), it is straightfor- 
ward to show that we can choose gi2 = gis = ^23 = 0, with gn, ^22, gss 7^ 0. 
This self-consistent choice leads to a diagonal mean-quadrupole tensor in a 
convenient laboratory frame of reference. We then write the quenched free 
energy 



3q = \ {ill + ^22 + ^33) j P W dXlnlJ dVl exp 



(3x J2 s'-'-mq,, 

fi=l,2,3 

(20) 



where 



/ dns"" (n) exp 



qtMti = Xp (A) dX- 



J^f^xs'^" m q,. 



J dfl exp 



J2/3XS>^>^ (Q) q,. 



(21) 
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Also, we remark that 

E^MM = 0, (22) 

which confirms the traceless property of the mean-quadrupole tensor. There- 
fore, we introduce the standard parametrization 

/giiO o\ I f v-s o\ 

q = g22 = - \, (23) 

V g33 / V 2s J 

so that (i) s 7^ and rj ^ in a. biaxial nematic phase, (ii) s ^ and rj = in 
a uniaxial nematic phase, and (iii) s = and 77 = in the disordered phase. 

The analysis of the quenched free energy depends on the choice of the 
distribution p (A). For example, we may choose 

p(A) = c5(A-l) + (l-c)(5(A + l), (24) 

which represents a sample with a number concentration c of prolate molecules 
(A = +1) and 1 — c of oblate molecules (A = —1), and which is convenient 
for comparisons with the annealed situation. If we adopt this form of p(A), 
it is straightforward to analyze equations ( l20l) and ( 12T]) . with the standard 
parametrization ( 123|) . and draw the phase diagram of figure (1). We indicate 
two uniaxial nematic phases, with s > 0, and with s < 0, separated 
by a first-order boundary (heavy dashed line) from the isotropic phase. The 
biaxial nematic region is limited by two critical lines that meet at the Landau 
multicritical point (c^ = 1/2 and = 1/(3l = 3/10). In the neighborhood of 
this Landau point, we confirm these results by the analysis of an expansion 
of the free energy in terms of the invariants of the tensor order parameter. 



2.2 Annealed disorder 

In the annealed approach, we write the canonical partition function 



{^i} {it.} 



2N 



E 



N 



1=1 



where the sum over the shape configurations {Aj} is restricted by the fixed 
value of the number density c of prolate molecules. 



N 



- = (2c - 1) . 



i=l 
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Figure 1: Phase diagram, in terms of the ration between temperature t/ti 
and the convcentration c, for the case of quenched disorder with a double- 
delta distribution. We indicate the nematic biaxial (S), two uniaxial nematic 
phases of opposite symmetry (A^+ and NJ) and the Landau multicritical point 
(L), located at = 3/10 and cl = 0.5. The heavy dashed line corresponds 
to first-order transitions. Heavy solid lines correspond to continuous tran- 
sitions (between biaxial an uniaxial nematic structures). We also indicate 
asymptotic results coming from an expansion of the free energy in the neigh- 
borhood ot the Landau multicritical point (thin solid lines). The stability 
curves for the isotropic and nematic phases are given by the dotted lines. 
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It is then convenient to introduce a chemical potential and to change to a 
grand canonical ensemble, 



N 



i=l 



+ 



2N 



E 



N 



i=l 



In analogy with the treatment of the quenched case, we use a polar 
parametrization for S'^", take advantage of the Gaussian identities to elimi- 
nate the squares, and write the asymptotic (A^ oo) result 



exp [-f3N(P] 



where is a grand potential per molecule. 



+ 2 (^n + ^22 + 933) - 1^ C 



with 



exp ( 



+ 



exp 



4"" 



and 



L {n, = S'' (Q) gn + S''' (Q) ^22 + S^' (Q) q,s 



(25) 



(26) 



(27) 



(28) 



which should be supplemented by the coupled equations for the minimization 
of (p with respect to qu, ^22 and ^33. Again, we see that the mean quadrupole 
tensor is traceless, so that we can use the standard parametrization of equa- 
tion 022]). 

The analysis of the free energy shows that the biaxial solution (s 7^ 0, 
7^ 0) is thermodynamically unstable (it is a kind of saddle-point instead of 
a minimum of 0). We then draw the phase diagram of figure (2), in terms 
of the thermodynamic field variables t = dimensionless temperature, 
and chemical potential /i. The dashed lines indicate first-order boundaries 
between the uniaxial nematic phases iV+ and A^_, and between the isotropic 
and each one of the nematic phases. The multicritical Landau point (at 
fi = and ti = 3/10) is just a simple triple point. We can also draw the 
phase diagram shown in figure (3), in terms of temperature and concentra- 
tion, which may be more interesting from the experimental point of view. 
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Figure 2: Phase diagram for the annealed case in terms of thermodynamic 
field variables (dimensionless temperature t = l/f3 and chemical potential 
fi). The triple point is located at = 3/10 and /i = 0. Dashed lines indicate 
first-order boundaries. 



The tie lines in the ordered region indicate the coexistence of two distinct 
uniaxial nematic phases. Again, we confirm these results by an analysis of 
an expansion of the free energy in terms of the invariants of the tensor order 
parameter. 

2.3 Two- temperature formalism 

In the two-temperature formalism, we introduce two heat baths, at differ- 
ent temperatures, associated with the relaxation times of the orientational 
(quicker) and disorder (slower) degrees of freedom. We now give a brief 
account of this formalism [17]. Given a configuration {Aj} of the slower dis- 
order variables, we can schematically write the probability of occurrence of 
a configuration {o"j} of the orientational variables, 

P \{\}) = ^ exp {{a,} , {A,})] , (29) 
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t/t. 




Figure 3: Phase diagram, in terms of temperature and concentration, in the 
annealed case. The tie hnes indicate the coexistence of two uniaxial ncmatic 
phases. The region of coexistence of the uniaxial nematic and the isotropic 
phases is too narrow to be represented in this graph. 



where T — 1//3 is the temperature of a heat bath, and 

Za = ({A,}) = [-/^^ ({^^} ' {^^})] ■ 

The time evolution of Aj is given by a Langevin equation. 



(30) 



dt 



-z (t) Aj - — + (t) , 



(31) 



where z (t) is a multiplier associated with the chemical potential, and 

{rji{t)rj,{t')) = 2rT^S,jS{r-t'), (32) 

where we have introduced the temperature Tx of a second heat bath. With 
the assumption of quick and slow time scales, it is reasonable to replace 
&H/d\i by its average value. 



m /dn\ _ meff 

dXi \dXi/^ dXi 
12 



(33) 



where 



■^e// = -Heff {{\^}) = -fc^Tlii ^^P [~ ^ ({a,} , {A,})] . (34) 

We then assume that the probabihty of a configuration {Aj} is given by the 
grand-canonical expression 



1 



exp [iSxfJ^Np - PxTieff] , 



(35) 



where 



E{/3,,/3,N,fi)= ! [rfA]i J]exp 

I {-} 



({c^.} , {A.}) + ^ ( EA. + iV 



(36) 

and the ratio n = T/Tx resembles the number of replicas in spin-glass prob- 
lems pi]. 

According to this two-temperature formalism, the orientational degrees of 
freedom of the Maier-Saupe model are allowed to reach thermal equilibrium 
at a temperature T and the disorder degrees of freedom thermalize at a 
temperature Tx 7^ T, with n = T/Tx- The problem is then reduced to 
the calculation of the grand-canonical partition function of n replicas of the 
original system [T6] . 



= Yl Yl { 5Z 



a=l 



2 



N 



N 



i=l 



2N 



\i=l 



In the thermodynamic limit, we write 

Etyjo ~ exp [-PN(t)two] , 

where 
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^two = "2/^^ + 2^ (^11 + ^22 + ^33) - ^ ^"^Ctwo 



(37) 
(38) 



with 



^two 



exp ( 



dnew[l3L{n,{q^^])] \ + 
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+ 



exp 



(39) 



The minimization of the grand potential iptwo leads to the equilibrium values 
922, and = —qu — 522- Note that we regain the annealed case for 
n = 1, and that the role of this parameter n will become clear in the next 
Section. 



3 Connections with the Landau-de Gennes 
expansion 

We have already chosen a standard order parameter, given by the traceless 
diagonal tensor fl23|) . We then introduce the second and third-order invari- 
ants, 

h = Tr [q2] = ql, + ql, + ql, = ^ {Ss' + rj') (40) 

and 

h = Tr [q3] = ql + ql, + ql, = {s' - rj') , (41) 

in terms of which it is usual to write the phenomenological Landau-de Gennes 
expansion for the free energy in the neighborhood of a transition, 

/ = /o + |/2 + f /3 + J ihf + jhh + I ihf + y (/2)' + ... (42) 

The Landau multicritical point is given hj A = B = 0. The stability con- 
ditions of the ordered phases in the neighborhood of the Landau point are 
discussed in terms of the signs of the coefficients of the higher-order terms. 
It has been shown that > is a necessary condition for the stability of a 
biaxial nematic phase in the vicinity of the Landau multicritical point p] |19]. 

In the present case, we remark that it is more convenient to adopt an 
alternative parametrization, in terms of two new variables r and ijj, such 
that 

gii = (^cos ip + a/3 sin , (43) 
Q'22 = (^cosV' - A/3sin'?/'j , (44) 

and 

Q'33 = rcosip. (45) 
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We then have qu + ^22 + 933 = 0, 

I2 = \r' (46) 

and 

/g = ^r=^ cos {"ill)) . (47) 
3.1 Quenched disorder 

In the neighborhood of the Landau muhicritical point, the quenched free 
energy, given by equation (l20l) . leads to the expansion 



5A4 „ 100A5 ^ ^ 2000A6 „ 1450A6 „ ^ 

+ ^3^2 + -^rihh + ^T^3 ^T^2 (48) 

42 (>?^ 231 (^A2) 7007 {>?^ 27027 (^As) 

where Qq is the free energy of the isotropic phase, 

A^= /" AV(A)rfA, (49) 



and we have kept terms up to fifth order only. The Landau multicritical 
point is given by t = BA^/IO and A'^ = (for an arbitrary distribution of 
shapes). 

Let us choose the double-delta distribution, given by equation (El]), which 
is particularly adequate for a comparison with the annealed case. The ex- 
pansion of the free energy in the neighborhood of the Landau multicritical 
point (/3 = /3l = 10/3 and c = cl = 1/2) is given by 

ln(47r) If 3/3 A , 10 , 

9a = -r^ + - 1 - — /2 2c - 1 /3+ 

/3 2 V lOy 21 ^ ' ^ 

1^2 , 2000^2 _ 245^^3 .5o^ 
42 ^ 7007 ^ 27027 ^' ^ ^ 

The positive sign of the coefficient of /| indicates that biaxial nematic phase 
is stable in the neighborhood of this Landau point. Sufficiently close to 
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the Landau point, we can show that the first-order transition between the 
uniaxial nematic and the disordered phase is asymptotically given by 



l^l.g(2c-l)^ (51) 

in agreement with numerical calculations (see figure 1). Also, we show that 
the critical lines separating the biaxial and the two uniaxial nematic phases 
are given by 

t 10 , a/3 /1001\^/^ , ,2/3 . ^ 

(6) ^ (tt^J (2c -1)'/% (52) 



tL 21' ' V1200, 

which also agrees with numerical calculations close to the Landau point (see 
figure 1). This phase diagram, with a stable biaxial nematic phase, is in 
qualitative agreement with previous results for a Maier-Saupe model with 
restricted orientations [13] [16]. Using the notation of the phenomenological 
Landau-de Gennes free energy, we remark that C > and = 0, as in the 
work of Allender and Longa [19]. However, new topologies may arise if we 
consider other forms of the distribution p (A). 



3.2 Annealed disorder 



In the annealed case, we use the grand potential 0, given by equation fl26|l . 
to locate the Landau multicritical point {(3l = 10/3 and fiL = 0), and write 
the expansion 

1 / 3f3\ ^ 50 ^ 
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5 ,2 2500 ,2 
42 2 27027 ^ 



63' 
1450 
27027 



'2' 



where 



1 ^1 



Svr cosh ( 



(53) 
(54) 



is the grand potential of the isotropic phase. The negative coefficient of J| 
shows that there is no stable biaxial nematic phase in the neighborhood of 
the Landau multicritical point. The line at = 0, below the temperature 
of the Landau point, is a first-order boundary between two distinct uniaxial 
nematic phase (with s > and s > 0). We show that the first-order lines 
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separating the uniaxial nematic from the isotropic phase are given by the 
asymptotic expression 

t 250 2 

in the immediate vicinity of the Landau point. The phase diagram in figure 
2 is in agreement with these asymptotic results. 

We can also calculate some asymptotic expressions in terms of tempera- 
ture and concentration, which are more convenient variables from the exper- 
imental point of view. For example, the region of coexistence of the uniaxial 
nematic phases in figure 3 is limited by the asymptotic border 

The asymptotic form of the first-order border between the uniaxial nematic 
and the isotropic phases is given by 

rr'^m{^-2)' '''' 

in full agreement with numerical calculations. Also, it should be remarked 
that the phase diagram in figure 3 is in qualitative agreement with previous 
calculations for a Maier-Saupe model with restricted orientations [TB] . 

3.3 Two-temperature formalism 

The same sort of calculations can be carried out for the grand potential 
in the two-temperature formalism. From equation (1381) . in the immediate 
neighborhood of the Landau point, we have the expansion 

1 3/3\ ^ 50 /i^ 

4>two = (PO,two + o 1 - 777 F^-'2 - -^-^3 + 

2 \ 10 J 63 n 

11/2 , ^500_^ /108 _ 143n\ _ i450_^^3 

42n2 ^ 27027n4 V 35 35 J ^ 27027 ^ ' 

where (pQ^two is the grand potential of the isotropic phase, and we keep terms 
up to fifth order. Of course, we recover the expansion for the annealed 
case with n = T/Tx = 1. The Landau multicritical point is still located at 
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— 10/3 and ji^ = 0, but the sign of the coefficient of /| depends on the 
parameter n. Indeed, there will be a stable biaxial nematic phase for 



< T7TT ~ 7T, (59) 

143 3 ^ ^ 

which indicates that a slight departure from complete annealing {n — 1) 
is already enough to give rise to a stable biaxial structure. Comparisons 
with previous results for a Maier-Saupe model with restricted orientations, 
n < 0.9, show that a somewhat larger difference between the temperatures is 
needed to stabilize the biaxial phase in the presence of additional direction 
fluctuations. 



4 Conclusions 

We have carried out exact statistical-mechanics calculations for a Maier- 
Saupe lattice model with the inclusion of extra disorder degrees of freedom 
to mimic a mixture of discs and cylinders. The closed-form solutions can be 
written as a Landau-de Gennes expansion for the free energy in the neighbor- 
hood of the transition, with explicit forms of model-dependent coefficients, 
which allows the use of several results from the literature. The stability of a 
biaxial nematic structure depends on the treatment of the disorder degrees 
of freedom. For quenched disorder, with a typical double-delta distribution 
of discs and cylinder, we obtain a global phase diagram, in terms of tempera- 
ture and concentration, with a Landau multicritical point, a biaxial and two 
uniaxial nematic phases. If the disorder degrees of freedom are allowed to 
reach thermal equilibrium, we show that the biaxial structure becomes un- 
stable. We then assume that orientation and disorder degrees of freedom are 
coupled to different heat reservoirs, with two different temperatures. In the 
two-temperature context, we show that a small temperature difference, which 
is equivalent to a small departure from thermalization, is already enough to 
stabilize a biaxial nematic structure. These results explain some disagree- 
ments in the earlier literature. Also, they qualitatively agree with a previous 
calculation for a similar Maier-Saupe model with restricted directional ori- 
entations. 
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